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Abstract 

We study the bead-spring model for a polymerized phantom membrane in the overdamped limit, 
which is the two-dimensional generalization of the well-known Rouse model for polymers. We derive 
the exact eigenmodes of the membrane dynamics (the "Rouse modes"). This allows us to obtain 
exact analytical expressions for virtually any equilibrium or dynamical quantity for the membrane. 
As examples we determine the radius of gyration, the mean square displacement of a tagged bead, 
and the autocorrelation function of the difference vector between two tagged beads. Interestingly, 
even in the presence of tensile forces of any magnitude the Rouse modes remain the exact eigenmodes 
for the membrane. With stronger forces the membrane becomes essentially flat, and does not get 
the opportunity to intersect itself; in such a situation our analysis provides a useful and exactly 
soluble approach to the dynamics for a realistic model flat membrane under tension. 

PACS numbers: 05.40.-a, 02.50.Ey, 36.20.-r, 82.35.Lr 
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I. INTRODUCTION 



Bead-spring models for linear polymers, wherein the polymer is modelled as a sequence 
of beads connected by springs, play a central role in the theory of polymer dynamics. For 
a linear bead-spring polymer chain, with position R„ of the n-th bead, n = 1 . . . N, the 

potential energy is given hj U = — (R-n. ~ R-n+i) ? in which k is the spring constant. 

2 „=i n 

Within the scope of bead-spring models, the Rouse model [1| holds a special status, as 
it provides a framework for linear polymers by which virtually any static and dynamical 
quantity of interest can be calculated analytically 2, In the Rouse model the dynamics of 
the beads is formulated in the overdamped limit. For solving the Rouse model one considers 
the so-called Rouse modes (rather than the equations of motion for the individual beads). 
Using the mode amplitude correlation functions, the quantities of interest for a phantom 
chain can be analytically tracked by reconstructing them from the modes {2, 3|. The Rouse 
model has been shown useful for self- avoiding linear polymers j4|, and to show that the 
dynamics of a tagged bead in a linear bead-spring model is described by the Generalized 
Langevin Equation (GLE) 

It is therefore a natural question to ask if the analytical solution of the Rouse model 
can be generalized to higher dimensions, e.g., making it useful for the dynamics of polymer- 
ized membranes. There is a considerable amount interest in the properties of polymerized 
membranes The interest stems not only from the point of view of a fundamental un- 
derstanding, but also because of their importance in biology and chemistry. Instances of 
polymerized membranes with a fixed connectivity can be found in 2D-cytoskeletons of cells 
SMlO] and graphite oxide sheets ll|-|l3[. 

Static properties of polymerized membranes as a two-dimensional extension of (one- 
dimensional) linear polymers came into fashion in the late 80s. At that time the driving 
question centered on the equilibrium properties of polymerized membranes, in particular, 
the scaling of the radius of gyration of a polymerized membrane with its lateral size, while 
the corresponding scaling behavior for phantom and self-avoiding linear polymers were al- 
ready well-known. Sophisticated renormalized group approaches were developed to this end 



14|-|l7]. accompanied by computer simulations 18|-|26|. A notable outcome of these studies 



is that the radius of gyration Rg of a membrane scales with its lateral size as i?^ ~ log A^ 
for phantom, and as Rg N for self-avoiding membranes. 
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In contrast to the equilibrium properties, the dynamical 
been studied with less intensity. Apart from scaling analyses 



properties of membranes have 



27H29| , the bulk of the research 



on the dynamics of polymerized membranes are heavily dominated by computer simulations 
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-|33|. leaving exact analytical results on the dynamics of membranes a relatively open area. 
In this paper we attempt to fill this void — we consider two-dimensional square-polymerized 
phantom membranes embedded in a three-dimensional space, both in the absence and the 
presence of tensile forces — and perform a Rouse mode analysis. We show that the Rouse 
modes are the exact eigenmodes of the membrane dynamics. Akin to the Rouse model for 
bead-spring linear polymers, our exercise allows us to exactly solve for the static and dynamic 
properties of phantom membranes; in this process also deriving the exact expression for the 



scaling of Rg as a function of the membrane's lateral size (which confirms the R 



ogA^ 



scaling), as well as show increased fluctuations at the edge of the membrane 3J, |35|. In 
particular, we draw the reader's attention to the exact solution under tensile forces, for which 
forces the membrane essentially encounters a fiat geometry. Unlike a phantom membrane 
in equilibrium, in such a situation the membrane becomes essentially fiat, does not get the 
opportunity to intersect itself, and therefore provides a highly useful and exactly soluble 
approach to the dynamics for a realistic model fiat membrane. We also note that the 
methods are generalizable to arbitrary internal and spatial dimensions. 

The structure of this paper is as follows. In Sec. [TTlwe describe the dynamical equations 
without forces, and their diagonalization by the Rouse modes. In Sec. IIIII we derive the 
scaling of the radius of gyration, the mean-square displacement of a tagged bead and the 
autocorrelation function of a vector connecting two beads. In Sec. IIVI we incorporate tensile 
forces in the dynamical equations and derive the new Rouse modes that diagonalize them. 
We finally end the paper with a discussion in Sec. |Vl 



II. DIAGONALIZATION OF EQUATIONS OF MOTION FOR THE POLYMER- 
IZED MEMBRANE AND THE ROUSE MODES 

We consider a rectangular polymerized membrane, for which N = Li x L2 beads are 
connected in a perpendicular toplogy, with Rn(t) denoting the spacial position at time t 
of the bead internally labeled as n. Naturally for a membrane it is convenient to take two 
numbers n = (ni, 712) for labeling the beads with Ui = 1, . . . , Lj. The potential energy for a 
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square-polymerized membrane is a simple extension of the Hamiltonian for a linear polymer 
in two internal dimensions 



^ Li— 1,L2 ^ Li,L2 — 1 

= ~ ^ ^ (-f^n.i,n2 ■R'?ii+l,n2) "I" ^ ^ ^ (R'ni,n2 -R-ni, 712+1; 
ni,n2=l ni,n2=l 



(1) 



for some spring constant /c. 

In the absence of externally applied forces, in the overdamped limit the dynamics of each 
bead of the membrane is described by 

dt C^Rn^"' 

where ( is the friction coefficient of the solvent and gn is the thermal force on the n-th bead. 
The thermal forces are uncorrelated between different beads, as well as in time, i.e.. 



(gm(t)-gn(t')) = ^W(^-0 

with Boltzmann constant ks and temperature T. 

For a linear polymer with positions R„(t) of beads n = 1, 
modes are given by [2|, |3| 



(3) 



N at time t the Rouse 



1 ^ 



cos 



n=l 



7i{n - l/2)p 



N 



R.n(^), 



(4) 



with p = 0, . . . , N — 1 and the inverse given by 

N-l 

R„(t) = Xo(t) + 2^cos 

p=i 



7i{n - l/2)p 
N 



X,(t) 



(5) 



The Rouse modes for the membrane are quite similar to those of the linear chain: since 
topologically the internal directions (connectivity) are orthogonal, the Rouse modes will be 
products of Rouse modes for a linear polymer. Below we introduce the following definitions 
for an elegant notation and proof of the independence of the Rouse modes. 
First we define 

2 if Pi = 1, . . . ,Lj - 1 
1 ifp^ = 



(3p 



/3p = (3r>A 



piyp2 



(6) 



and 



fp^{ni) = COS 



npi^rii - 1/2) 



/p(n) = /pi(rii)/p2(n2) 



(7) 
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The orthogonality relation at the basis of the proof for the Rouse modes for a linear polymer 
chain can then be generalized for the membrane 



(8) 



where the summation is taken over all allowed values of n. The Rouse modes amplitudes 
and the corresponding inverse are then given by 



RnW = J]/3p/p(n)Xp(t). 



(9) 
(10) 



The position of the center-of-mass Rcm(^) at time t simply equals Xo(t). 

Next, we note that when the equations of motion from Eq. ([2]) for the beads are expressed 
one gets a term for each internal dimension equal to that of a term for a bead in a linear 
polymer. Taking the time derivative in Eq. ([9]) and plugging in the equations of motion for 
the beads and then the inverse we get 

dXp(t) 



dt 



-apXp(t) + Gp(t) 



with ttr 









TTP2 




^sin^ 


2Li 


+ sin^ 


2L2_ 


) 



(11) 

(12) 



i.e., the Rouse modes diagonalizes the equations of motion. This set of linearized differential 
equations can be solved just as for a linear polymer by using the correlation function for the 
thermal forces given by Eq. By using the orthogonality relation in Eq. ([S]) 



(Gp(t)-G,(0) = ^W(^-0, 
so that the following two-point correlation functions for the eigenmodes can be derived 

l2\ _ /FTD /^N T3 ^nM2v QksT 



(13) 



Xoo(t) = ([Xo(t) - Xo(0)]') = ([Rcm(t) - Rc 



-t 



Xpq(t) = (Xp(t) ■ Xq(0)) = — ^- exp [-apt] 6^^ with p 0. 



(14) 
(15) 



CN(3p ap 

Just like the case of a linear polymer, any correlation function of interest for the membrane 
can be derived using Eqs. f |TllfT3l) . We address a few of them in the following section. 
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III. A FEW PROPERTIES OF INTEREST FOR A MEMBRANE 



All throughout this section, for the sake of simplicity, we restrict ourselves to square 
membranes with Li = L2 = L, and thus N = L x L. 



A. Scaling of the radius of gyration Rg 

The radius of gyration squared for the membrane is defined by 

n 

= 1^ E E E(^p(^) ■ Xq(t))/3p/3Jp(n)/,(n), (16) 

n pt^Oqt^O 

in which Eq. (ITO!) is used for the second line. Having plugged in Eqs. ([8]) and (fTSj) the 
squared radius of gyration reduces to 

3kBT 1 

The summation (ITTj) . using the definition of Op (fT2|) . can be split into two parts: 

3k l\ knT 

The first term in Eq. (1181) can be evaluated exactly: X = (A^ — 1)/A^. For the second term we 
take the long polymer limit ^ 1 so that a-p can be expanded and the sum approximated 
by an integral with Xj = pi/ L where the integral is over a unit square in the positive quadrant 
excluding the small area near the origin. The second approximation we make is replacing 
the area of integration by that of the positive quadrant of a unit circle excluding the small 



Rg 




L-1 L-1 L-1 

Tip' 



Pl = lp2 = l 



area near the origin which can be solved analytically 

c/x 3 dr 



^ = ^ / ^ = ^ / — = ^ logiV. (19) 
The radius of gyration squared for a square membrane in the long polymer limit thus becomes 

Ji^^.tff^,^^,,^). (20) 



^ k \ N ' An J 
A comparison between the two terms in Eq. f|20|] shows that the second one is bigger than 
the first one for > 61. In Fig. [T]we present a comparison between the analytical result 
(120]) and the exact evaluation of Eq. ( IT7|) for L up to L = 10^. Our result (I2U1) confirms the 



log scaling of i?^ obtained earlier by field-theory methods 14|-|l7|. 
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FIG. 1: The radius of gyration squared i?^ for a square membrane as a function oi N = L x L 
beads, with ksT/k set to unity. The black dots represents the sum in Eq. (|17p evaluated exactly 
for L up to L = 10^, while the continuous red line represents the function (iV — 1) /A^ + 3/(47r) log N. 

B. Mean-square displacement of a tagged bead 

Let us define ARn(t) = Rn(^) ~ H-nlO) for the n-th bead. We use Eq. (fTO|) to write 

(AR„(t)2) =([Xo(t)-Xo(0)]2) 

+ 2 E E /3p/9q/p(n)/q(n) {Xpq(O) - Xpq(t)} . (21) 

Using Eqs. (I14til5p this simplifies to 

At short times kt/( ^1/8 the exponential in the second term can be expanded, and the 
sum is exactly evaluated. In comparison to the second term the first one can be neglected 
to obtain 

(AR„(t)^) = E ^Mnr = (23) 

which agrees with the mean-square displacement of a free bead, and confirms the wisdom 
that at short times the beads do not feel the connectivity. 
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At very long times the exponentials in Eq. (122|) reduce to zero, the second term 
becomes a fixed number independent of t [but dependent on the location of the bead: 
the sum is largest when the bead is located at the corner of the membrane, and equals 
QkeT (logiV/vr - 7r/4 + 2/3) /% i.e., for kt/C > iVlogiV/7r the first term of Eq. 1^2) domi- 
nates, and the motion of the tagged bead becomes simply diffusive, with the same diffusion 
coefficient as that of the center-of-mass. 

At intermediate times the mean-square displacement of the tagged bead will depend on 
its location on the membrane. We work out three difference cases when the tagged bead 
is (a) the central bead of the membrane, (b) located in the middle of an edge, and (c) a 
corner bead. The enhanced mobility of edge and corner beads was already reported for 
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35| 



self-avoiding membranes 

(a) For the mean-square displacement of a bead at the center of the membrane (ARm(t)^) 
the sum in Eq. (122!) will only be over even values of pi. In the long polymer limit the sum 
can be converted to an integral with Xi = Pi/L, such that 



(AR„,(t)2) = ^ f'dx, /'dx2— {l-e-"pn = ^ f'dt' C dx^ C dx^e 
C io io «P C io Jo Jo 



(24) 

We proceed with Eq. (12^ by splitting the exponentials, resulting in two identical integrals 
which can be evaluated in terms of the modified Bessel function of the first kind Ia{z) and 
making a change of variables. 



(AR„,(t)2) 



dt' 



dxexp 



Akt' 



sm 



C Jo \Jo I C 

To characterize the moderate time behavior of Eq. 



vr ■ 

—X 

2 J 



GknT 
k 



dt'e-^''ll{2t'). 

(25) 



we look at the asymptotic expansion 



Ia{z) = ex]){z) / ^/27TZ and evaluated the integral so that 



kt 

X 



(26) 



4.04. 



(27) 



where 

C^= lim -log[a;]+47r / dte-^^I^{2t) 

(b-c) The same can be done for the mean-square displacement for a bead located in the 
middle of an edge (ARe(t)^) and corner (ARc(t)^). The corresponding calculations are 
similar to those in (a), and they yield 

^kf 



(ARe(t)^) = ^ ( log 



c 



(28) 



where 



and 



where 



Ce = hm - log[a;] + 2ti / dte"^* (/o(2t) + h{2t)) Io{2t) ^ 2.47, 



(ARe(t)^) = ^ ( log 



C, = hm - log[x] + vr / rfte"^* (/o(2t) + /i(2t))^ f« 1.47. 



(29) 
(30) 
(31) 
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FIG. 2: The mean-squared displacement (ARn(t)^) for beads in a square membrane as a function 
of scaled time kt/Q. The sum in Eq. (|22p was exactly evaluated for ksT/k = 1 for beads at the 
center of the membrane, at the middle of an edge, and at a corner for N = 301^ and the result is 
represented by the circular (red), triangular (green), and square (blue) data points respectively. The 
cyan diamonds represent the exact evaluation for a bead with internal position n = (151, 5) to show 
the transition that the mean-squared displacement for a bead makes from one region to another. 
For very short times kt/C, ^ 1/8 the mean-squared displacement for all beads (ARn(t)^) ~ t thus 
behaving like that of a free bead as shown in Eq. (|23p which is represented by the dashed line in (a). 
For very long times the mean-squared displacement for all beads (ARn(t)^) ~ t/N behaves like that 
of a mean-squared displacement of the center-of-mass which is represented by the dotted line in (a). 
In the intermediate time regime the mean-square displacement for beads (ARn(t)^) ~ log[A:t/C] 
as in agreement with Eqs. (j26p . (j28p and (p0|) — the corresponding approximations are valid for 
kt/Q Nl-n"^ — and are represented by the red, green, and blue solid lines in (b) respectively. 
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Note that all constants C"s can be evaluated up to arbitrary precision, and that that the 
approximations ([25]), ([27]) and ([29]) are valid for kt/C, < N/tx'^. In Fig. |2] we compare the 
exact evaluations from Eq. ([22]) and the corresponding approximations ([23]) . ( 126]) . ([28]) and 



C. Autocorrelation function of a vector connecting two beads 

In the most general case the vector rjnn(^) = R-ml^) R-nl^) connects two beads with 
internal coordinates m and n at time t. The autocorrelation function of this vector -Dmn(^) = 
(rmn(i) ■ rmn(O)) cau bc expressed in terms of the eigenmodes 



p^O "P 



Like in the case for the mean-square displacement of a tagged bead, the behavior of -Dmn(^) 
depends on time and the internal positions of the beads, hence we only briefly outline the 
calculation procedure and state the results. We focus on two extreme cases (a) D(.{t) where 
the vector connecting the two beads located at the opposite corners of the membrane, and 
(b) -Da(i) where the vector connects two beads at position n and n + a somewhere in the 
middle of the membrane for a/L = |a|/L ^ 1. We will generalize these results qualitatively 
for the cases when the locations of the beads and the distances between them are arbitrary. 

Short times: For times kt/C^ <^ 1/8 the exponent is roughly 1 and so the functions 
are constant. In case two beads are close neighbors the sum can be evaluated -Da(^) = 
kBT{3/2 + loga)/k, which is exact for a = 1, and a rough estimate Dc{t) = GksT log N/ (irk) 
follows from a calculation similar to that performed for the radius of gyration. 

Intermediate times: For a/L ^ 1 the sum for Dg^{t) can be expanded in a and evaluated 

as 

Further, it can be shown that the summation for Dc{t) equals that of — (ARc(t)^) up to a 
constant, such that 

. (34) 



ah, 

D,{t) = DM - ( log 



kt 
J. 

Long times: For a/L <^ 1 and kt/( ^ N/ti"^ only the lowest mode contributes to the 
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FIG. 3: The scaled autocorrelation function -Dmn(i)/-Dnin(0) of the vector connecting beads with 
internal position m and n in a square membrane as a function of scaled time kt/C,. The sum in 
Eq. (|32p was exactly evaluated with ksT/k = 1 for two beads at positions [(251, 251), (250, 251)], 
[(251, 254), (258, 251)], [(1,1), (501,501)], and [(251, 281), (281, 251)] in the membrane with N = 501^, 
with the results shown by the circular (red), triangular (green), square (blue), and diamond (cyan) 
data points respectively. Note that the first two pairs of beads are very close to each other in the 
middle of the membrane where the third pair are two beads in opposing corners of the membrane; 
i.e., the first two cases correspond to Dg^{t)/ Dg^{0), while the third one corresponds to Dc{t)/ Dc{0). 
The function D{t) ~ log[kt/(^], as in agreement with Eq. ()34p corresponding to the blue solid line 
in (b), for a duration that is larger when the beads are connected by a larger vector a. For very 
small vectors a the behavior will shift to D{t) ~ 1/t as derived in Eq. (|33|) corresponding to the 
red and green solid lines in (a) before kt/C ^ N/tt'^. The function D{t) for the last pair shows a 
transition from logarithmic behavior to D{t) ~ 1/t before the diffusive time regime. 



summation and so expanding the correlation function in Eq. (132|) yields 

TT^k 



while 



QknTa'^ 



D.(t) = -^exp 



(35) 



(36) 



It is interesting to note the differences in behavior for -Da(^) and Dc{t) at intermediate 
times. The reason behind this difference is as follows. For -Da(^) the beads are close and 
as a result they quickly become 'aware' of each other's presence. For Dc{t) on the other 
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hand, the beads do not become 'aware' of each other's presence almost until kt/C, ~ N/n"^. 
Based on this observation we expect that when the beads are neither very close nor very far, 
there will be first a logarithmic decay (IMI) for -Dmn(^); followed by a 1/t power-law decay 
(155]) before the terminal exponential decay fl35ti36p sets in. The above results are verified 
by comparing the approximations to the exact evaluation of Eq. (15^ for different pairs of 
beads in Fig. |3l 



IV. DYNAMICAL EIGENMODES IN THE CASE OF EXTERNAL TENSILE 
FORCES 

In this section we concentrate on presenting the exact eigenmodes when the membrane 
is stretched in two perpendicular directions by forces Fi and F2 applied at the edges of the 
membrane. 

Adding tensile forces to the system ensures that each bead has its own mean position 
around which it fluctuates due to thermal forces. For large enough tensile forces the beads 
are far enough away from each other, and self-intersection of the membrane is avoided. Such 
a situation therefore mimics the behavior of a realistic flexible membrane under tension. 
The system is still exactly solvable by introducing the following term to the Hamiltonian 

Li-1,L2 ii,L2-l 

Uf= Fi-(R 

rii,n2 -Rni+l,n2 

(37) 

ni,n2=l ni,n2=l 

making the exercise of this section useful for practical purposes. 

The mode amplitudes and their inverses are once again defined by Eqs. (I WTU]) . but the 
following term Hp will be added to the right hand side of the differential equation in Eq. 
(fTT]) in order to solve for the dynamics of the mode amplitudes: 

= ^Xl^P^'^)[^l^'^"i^i"^"il)+^2(5n2L2-5„2l)] 



F F 

^P2o{/pi(^i) - /pi(l)}7^ + ^pio{/p2(^2) - /p2(l)}7^ 



(38) 



Note that spatial symmetry for the modes is broken, and that Eq. ( |T5l) is replaced by similar 
equations but also depending on the spatial components i and j of the vectors. 

Xp,^(t) = (Xp,(t)Xq,(o)) = + [-"p^] ^pq^^^- p ^ 0. (39) 

"p C^Pp "p 
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The assumption Fi _L F2 is necessary in order for the modes to be the exact dynamical 
eigenmodes. If the tensile forces are not orthogonal, then it will result in correlations between 
some of the modes (and correspondingly a parallelogram- like structure of the membrane). 




FIG. 4: A randomly generated membrane consisting of 25 x 25 beads in a three dimensional space, 
with orthonormal tensile forces, ksT = 1/2, and k = ( = 1. The black lines connects neighboring 
beads. Due to tensile forces the membrane gets a rectangular lattice like structure where thermal 
forces push the beads out of their mean positions. 



It is of course interesting to investigate the transverse fluctuations of the membrane under 
lateral tension. Two observables we consider are the effective thickness and the additional 
surface area under thermal undulations. A quantity that can be used as a measure for the 
thickness of the membrane, is the standard deviation of the relative height of the bead in 
the middle of the membrane with respect to its mean position. In the limit of a continuous 
very large membrane, this thickness is given by 



Ak ■ 



(40) 



We verified numerically that the fluctuations increase for beads near the edge of the mem- 



brane, in line with self-avoiding membranes 



34 



35|. Note that this thickness is unaffected 
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by the strength of the tensile forces. 

A much more involved calculation is needed for the area of the membrane. Consider a very 
large square membrane with N = L"^ beads and perpendicular tensile forces Fi = F2 = F 
of equal strength. We consider the case where the temperature is low enough, or the tensile 
forces are strong enough, such that the fluctuations of the beads around their equilibrium 
are small compared to the distances between two neighboring beads. The total area of 
the membrane can be written as the sum of areas of all triangles between the three beads 
with indices {i + 1, j) and (i, j + 1) plus the sum of areas of all triangles with indices 

(z + 1, j), + 1) and {i + 1, j + 1); The expectation value for this area, in the continuum 
limit, is then given by 



- 1)2 kBT{2L^ - L) 

Equations f l4mHT]) collectively imply that under the application of tensile forces, irre- 
spective of the force magnitude, the out-of-plane degrees of freedom decouple from in plane 
degrees of freedom. Moreover at strong stretching F ^ ksT/k the main contribution to 
the area of the membrane comes from the stretching of the individual springs. As already 
mentioned earlier, at strong stretching there is no difference between a phantom and a 
self-avoiding membrane, and therefore we expect the same for the increase in area for a 
self-avoiding membrane under strong stretching forces. 

The eigenmodes also give access to dynamical information. For instance, to obtain the 
mean vector connecting two neighboring beads in the membrane ARi(t) = R„i-|_i^„2(^) ~ 
^ni,n2{'t) with rii = 1, . . . , Li — 1 and 71,2 = 1, ... , L2, we use (Xp(t)) = Hp/ctp, which can 
be proved by solving the differential equation in a similar fashion as for Eq. ( l39l) . such that 

(ARi(t)) = V/3p^^^[/p(ni + 1,^2) - /p(ni,n2)], (42) 

yielding ARi(t) = Fi/k. The same is true for neighbors in the other internal direction so 
that for orthogonal tensile forces, on average, the membrane obtains is a fiat rectangular 
structure. (The above results are obtained by using 

L-l 



^sin 1 

p=i 



'TVp' 




'irnp' 




'np' 




sin 




cos 




[2l\ 


[ L \ 


[2l\ 



I forn = l,...,L-l, (43) 



where the summation is over odd values of p). Further, the random thermal forces are 
gaussian distributed with a variation given by Eq. ([3]) such that (Xpj(O)) and (Xp(0)) can 
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be used to derive the first two moments for every mode amplitude. We use tliese results to 
generate a random typical configuration for a membrane as shown in Figure |H 

V. DISCUSSION 

The Rouse model has proved to be a cornerstone for polymer dynamics. We generalized 
this model with two internal coordinates, for which the entity becomes a phantom mem- 
brane. We have shown that a set of eigenmodes can be used to solve the membrane dynamics 
in the overdamped limit, and have demonstrated that many interesting properties can be 
analytically derived for the membrane using certain relations for these eigenmodes. Fur- 
thermore, we have shown that adding large enough tensile forces to the system mimics the 
behavior of a realistic flexible membrane under tension: the eigenmodes can still be exactly 
solved analytically, making the analysis useful for practical purposes. 

We also note, although the following is not relevant to study the properties of a phantom 
membrane, that the entire analysis presented in this paper is easily generalized for any 
number d > 1 of internal coordinates. For every internal coordinate the Rouse modes and 
inverse get a factor fp.{mi), /3p is a product over d factors, and Op becomes a summation 
over d terms. By generalizing the potential energy so that all neighbors are connected, Eqs. 
flSl lTS]) remain valid in their current forms. It can also be shown that adding a tensile force for 
each internal coordinate still makes the system analytically solvable. Similarly, the analysis 
can also be extended to incorporate periodic boundaries (resulting in a torus or cylinder 
geometry, as appropriate, when the internal dimension d = 2). 
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